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Abstract. Synchrotron diffuse radiation (SDR) emission is one of the major Galactic components, in the 100 
MHz up to 100 GHz frequency range. Its spectrum and sky map provide valuable measure of the galactic cosmic 
ray electrons (GCRE) in the relevant energy range, as well as of the strength and structure of the Galactic 
magnetic fields (GMF), both regular and random ones. This emission is an astrophysical sky foreground for 
the study of the Cosmic Microwave Background (CMB), and the extragalactic microwave measurements, and 
it needs to be modelled as better as possible. In this regard, in order to get an accurate description of the SDR 
in the Galaxy, we use - for the first time in this context - 3-dimensional GCRE models obtained by running the 
Dragon code. This allows us to account for a realistic spiral arm pattern of the source distribution, demanded 
to get a self-consistent treatment of all relevant energy losses influencing the final synchrotron spectrum. 


1 Introduction 

Defiection of ultra-high energy cosmic rays (UHECR), rota¬ 
tion measure, synchrotron radiation, and polarized dust are 
just a small sample of different methods of observation able 
to investigate the galactic magnetized interstellar medium 
(ISM). In this frame, cosmic rays (CRs) are, doubtless, a 
unique probe of the ISM properties. Thanks to a set of suc¬ 
cessful experiments such as Fermi-LAT, PAMELA, AMS- 
02, the last few years have witnessed an incredible progress 
in the knowledge of electron and positron Galactic CRs, over 
a wide range of energy, from O(TeV) down to tens of MeV. 
Unfortunately, solar modulation complicates matters, since 
the CR spectra observed on Earth are - for < 20 GeV - 
completely reshaped with respect to their local interstellar 
spectra (LIS). 

Relativistic cosmic ray electrons and positrons (CRE), spi¬ 
ralling around the interstellar magnetic field lines, are at the 
origin of the radio diffuse emission from the Milky Way. For 
magnetic field intensity of 0(/iG), like in the case of our 
Galaxy, and for electrons/positrons of [GeV^TeV] energies, 
the synchrotron emission falls in the [MHz-^GHz] range. In¬ 


deed, the SDR offers valuable complementary checks of the 
low energy spectrum, and in general of the spatial distribu¬ 
tion of CRs in the Galaxy. Hence, a parallel study of radio 
emission, together with CR measurements, can put better 
constraints on all the interstellar medium (ISM) components 
involved ( [Strong et al.[ [201 Ij ). The interpretation of those 
measurements requires a proper modelling of injection, prop¬ 
agation and losses in the Galaxy. 

Moreover, the presence in the [20-^200] GHz range of sev¬ 
eral astrophysical sky signal components - with similar inten¬ 
sities and some spatial correlation - makes the extraction of 
the CMB a complex task. In order to achieve sufficient ac¬ 
curacy on the cosmological signal the component separation 
needs to take advantage of the knowledge on the properties 
of diffuse Galactic emission. 

We plan to accomplish the aforementioned study by run¬ 
ning the Dragon code in its 3-dimensional version. Indeed, 
this is well suited to model the CRE propagation, when ac¬ 
counting for a realistic spiral arm distribution of astrophysi¬ 
cal sources, gas distributions, magnetic fields models and dif¬ 
ferent position-dependent models for diffusion in the parallel 
and perpendicular directions with respect to the GMF. 
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Table 1 . The CRE models considered in the present analysis. The 
reported values of 7inj(e“) refer to energies below/above 4 GeV. 


Model 

(5 

UA(kms ^) 

p 

7inj (e ) 

KRA 

0.5 

15 

-0.4 

1.6/2.5 

KOL 

0.33 

35 

1.0 

1.6/2.5 


2 Objectives and method 


In the present Section, we outline the guidelines of the multi¬ 
wavelength analysis we have performed, in order to model 
the CRE spectra consistently with the diffuse synchrotron 
emission of the Galaxy. One of our main aims has been: 


1. To explore the physical properties - injection and propa¬ 
gation - of the local interstellar spectrum (LIS) of CRE, 
in the realm of low energies (< 7 GeV), by combin¬ 
ing the latest e~ and e+ measurements with the diffuse 
Galaxy radio emission, between 10 MHz and few GHz. 
Below that energy, we modelled the LIS of e+ + e~ on 
the basis of the observed synchrotron spectrum of the 
Galaxy, which is unaffected by propagation in the he¬ 
liosphere (see also e.g., [Orlando and Strong] ( |20T3] )); 


2. In parallel to that, the current study has pushed us to 
give an important constrain on the vertical scale height 
of the diffusion region in the Galaxy, by looking simul¬ 
taneously at the radio spectrum, the latitude profile of 
the synchrotron emission, and the positron fraction at 
energies below ^ 5 GeV. 


The structure of the GME is still not well understood. Gen¬ 
erally, a realistic, and accurate description of the synchrotron 
emission, as well as of its angular distribution, requires to 
take in account two main components for the GME: the reg¬ 
ular and turbulent ones. Regarding the ordered one, here we 
rely on a recent model, based on a wide and updated com¬ 
pilation of Earaday rotation measurements ( jPshirkov et al. 
|2011| ). It consists of two different components: a disc field, 
with a magnitude in the vicinity of the solar system, Bq, 
taken to be 2 /iG, and a toroidal halo field, with a thick¬ 
ness of ^ [0.2 ^ 0.4] kpc, and extended - above, and below 
- out of the galactic plane (GP) for [1 ^ 2] kpcQ A new, 
and much-improved model for the regular GME has been re¬ 
cently brought to the attention of the community ( jJansson 
[and Earrar||2012j ). Besides a disc field and an extended halo 
field, the peculiarity of this GME model is its X-shape in the 
r — z directions. The main implications based on such a or¬ 
dered GME model will be addressed in our future work. 

^We have included the regular field in order to make our model 
compatible with the current information. However, in our analysis, 
we have checked that only the halo component (Bhaio = 4 jiG) for 
the regular GME plays, albeit marginal, a role. 



Figure 1. Normalized electron density profile along the radial di¬ 
rection. The black line corresponds to a 2-dimensional, smooth CR. 
The red one corresponds to a 3-dimensional spiral arm pattern of 
sources. 


Instead, for what concerns the random component - actu¬ 
ally, the main responsible for the diffusion of charged par¬ 
ticles in the ISM - there is still a poor knowledge about its 
geometrical structure. As in [Han et alT] ( |2004| ), we have as¬ 
sumed it to fill a thick disk, modelled with an exponential 
vertical profile, and an effective scale-height accordingly 
to the equation: 


-^ran(^) — -^ran(0)6^P( ^h)' 


( 1 ) 


Eor the aforementioned purposes, we first run DRAGON, a 
new numerical package aiming to solve the diffusion equa¬ 


tion of CRs in the Galaxy environment (Evoli et al. 2008 
[Gaggero et al.||201^ . 

Conversely to a too simplified assumption of a diffusion 
spatially uniform in the thick disc, here in this contribution 
we account for a possible spatial dependence of the diffusion 
coefficient. 


D(p,R,z)=DoP^f(z) ^ 
\po 


( 2 ) 


p being the rigidity of the particle, P the particle speed in 
units of speed of light c, and f{z) indicates the spatial de¬ 
pendence of the diffusion coefficient. As predicted by the 
quasi-linear theory (QLT), that should be related to the fluc¬ 
tuating magnetic field, and hence as D{z)~^ (x B^q,^{z) oc 
e-xp{-z/zh). 

In this paper, for sake of simplicity, only two represen¬ 
tative classes of propagation regimes have been taken in 
consideration: the KRA (Kraichnan), and the KOL (Kol¬ 
mogorov). The main parameter characterizing those two 
models are reported in Table 

Eor each of them, we varied the scale-height of the diffu¬ 
sive halo in the range = [1 ^ 16] kpc, and the main dif¬ 
fusive parameters were determined in order to to minimize 
the combined against the boron-to-carbon ratio and the 
proton observed spectra. 
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Finally, at high energies (> 7 GeV), unlike our previous 
results presented in |Di Bernardo et aL] ( |2013| ), here we fix the 
spectral index, and the normalization of the injection spec¬ 
trum of the primary electrons and of the extra-component 
by tuning our models against the new data, as recently re¬ 
leased by PAMELA and AMS-02 collaborations, respec¬ 
tively, rather than on e“ + e+ spectrum measured by Fermi- 
LAT ( [Gaggero et al.[ |2014| ). With reference to the Figure 
we want to make it clear that the assumption of a sim¬ 
ple power-law (PL) distribution for energetic electrons and 
positrons, whose sources are smoothly distributed in the en¬ 
tire Galactic disc, leads to large overestimation of their en¬ 
ergy densities in comparison with the values deduced when a 
3-dimensional spiral arm distributions are used. In our opin¬ 
ion, the impact of this more realistic modelling of the particle 
distribution on the final synchrotron spectral maps is a cru¬ 
cial issue (see also the comparison between the Figures |^& 


3 The synchrotron emission of the Gaiaxy 


It is known that the synchrotron intensity depends on the spa¬ 
tial, and energetic distribution of CRE density, Ue, and the 
strength of the magnetic field perpendicular to the line 
of sight (LOS) to the observer. Once the CRE densities are 
computed by DRAGON - at all points of the computational 
grid - we use GammasKY to get the emissivities (i.e. power 
per unit volume per unit frequency per unit solid angle), for 
the regular and random fields, according to the standard for¬ 
malism ( |Longair[|2011| ), after having integrated over the par¬ 
ticle energy. 

The emissivity (in erg s“^ Hz“^) - in an uniform magnetic 
field - is partially linearly polarizec0 and usually subdivided 
in two components, respectively parallel and perpen¬ 
dicular to B± = B{r)sma, with a the angle between the 
direction of the magnetic field and the LOS. The polarization 
formulation will be used in our future work. Here, we show 
results based only on the total intensity, given by 

e{u,r) = VS-^B_i(r)F{x); (3) 


where x = vli/c^, being = (3/47r)(e/mc)5x7^ the 
critical synchrotron frequency, 7 the particle (electron or 
positron) Lorentz factor, and F{x) is defined in terms of 
Bessel functions. In the case of a randomly oriented magnetic 
field, the expected isotropic emissivity is computed accord¬ 
ing to|Ghisellini et al.|(ri988|). 



Figure 2. The average synchrotron spectra, for Zh — ^ kpc. We 
show the spectra obtained with (solid lines) and without (dashed 
lines) the spectral break in the e~ source spectra. Dotted lines 
are the contribution of secondary e~ source spectra. The contribu¬ 
tion of the regular GMF is shown as the dot-dashed line. The The 
normalization required for the random component field strength is 
Bran(O) = 7.6 /XG. 



Figure 3. Normalization of random GMF the vertical scale 
height. The 3(5) cr allowed by RM are represented in grey (light 
grey) bands. The red squares are the values used in our KRA mod¬ 
els in order to reproduce the observed spectrum at 408 MHz. 


Galactic synchrotron spectrum. We take care of correctly re¬ 
producing the observed 408 MHz radio synchrotron radia¬ 
tion as in |Haslam et al.| ( |1982| ), by tuning - time to time - 
the normalization value for the turbulent component of the 
GMF (see the Eq. [T]). We get the sky maps in GAMMAS KY 
by integrating the Galactic emissivity along the LOS, 



l.o.s 


(4) 


3.1 The total synchrotron intensity 

Given a GMF model, and for the representative CRE den¬ 
sity models aforementioned, the next step is to compute the 

^for a monochromatic and isotropic distribution of CRE. 


where e(z/,r) is the total emissivity given by the Eqs. In 
Figure we refer to the observed brightness temperature (in 
K), defined as ( [Rybicki and Lightman[|1986| ) 


T(z.) 


c^I{iy) 


(5) 
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The sky maps are subdivided into equal ar ea pixels fol- 


lowing the HEALPi)|^pixelization scheme of Gorski et al. 


(2005). We average the flux over the sky regions 40° < / < 
340°, 10° <b< 45°, -45° < 6 < -10°, where I and b are 
Galactic longitude and latitude respectively. We restrict the 
analysis to the regions out of the Galactic plane - but avoiding 
the polar regions - being the contamination from point-like 
and local extended sources expected to be the smallest. In 
addition to that, we avoid absorption effect at radio frequen¬ 
cies, and free-free emission at higher frequencies ( |Rybicki| 
and Lightman[ 1 19 8^ . Therefore, the observed Galactic dif¬ 


fuse emission in the radio band is, almost entirely, due to the 
synchrotron radiation of CRE moving errantly in the GMF. 

From tens of MHz to 23 GHz (and up to 94 GHz for 
Wmap), in such sky region we directly compare our sim¬ 
ulated models with the synchrotron spectra measured by a 
wide set of radio surveys at 22, 45, 408, 1420, 2326 MHz as 
well as Wmap satellite data at 23, 33, 41, 61 and 94 GHz. 

With reference to the Figure it is immediate to realize 
that the radio data (< 0(100) MHz) are clearly incompatible 
with a single PL electron spectrum, suitable to fit the CRE 
data. Rather, we find that introducing - below a few GeV - 
either a break or an exponential infrared (IR) cut-off in the 
population of primary e“, that helps us in providing a very 
good description of the radio data. 

As immediate consequence of that, the total electron flux 
(at < 4 GeV), and hence the radio spectrum below 100 
MHz, are dominated by secondary particles, which are pro¬ 
duced in nuclear collision with the nuclei of the gas present 
in the ISM, offering thus a direct probe of the interstellar pro¬ 
ton spectrum. In this regard, we found that once the low en¬ 
ergy e~ source spectrum is tuned to reproduce the observed 
e+ spectrum, only models featuring low re-acceleration can 
reproduce the observed e+ spectrum and fraction, in total 
agreement with what found in [Strong et al. | ( [MTT] ) ; [Orlando 
and Stron^ ( [2013| ). 


3.2 The magnetic halo height 

The vertical - perpendicular to the Galactic plane - size of the 
CR diffusion region represents, undoubtedly, a cornerstone 
in modem Astro-particle physics. The accurate knowledge 
of it goes beyond the target of conventional CR astrophysics; 
it is also worth for Dark Matter (DM) indirect search, since 
that the local flux of DM decay, and annihilation products are 
expected to depend significantly on such physical length. 

So far, the diffusive vertical boundary has been con¬ 
strained purely on the basis of CR radionuclide, 
ratio most commonly. However, this method is marked by 
several flawless, due to the severe uncertainties connected to 
local distribution of sources, gas, and especially by the solar 
modulation. 


^http://healpix.jpl.nasa.gov/ 


Sky region: 40° < I < 340°, ® v = 408 MHz 



Radio data e 40° < I < 100° 
KRA Model, zt = 2 kpc 
KRA Model, zt = 4 kpc 
KRA Model, zt = 8 kpc 
KRA Model, zt = 16 kpc 

-20 0 20 

Galactic latitude [deg] 


Figure 4. The latitude profile for the synchrotron emission at 408 
MHz, at different magnetic halo height. 


To the contrary, the synchrotron emissivity of the Galaxy 
offers a much more genuine probe of the scale height Zh- In 
that regards, we may notice that, when a realistic vertical dis¬ 
tribution is adopted for the radiation interstellar field (ISRF) 
and for the GMF, energy losses in the O(GeV) energy range 
- hence in the radio energy band - do not affect significantly 
the CRE vertical distribution, determined predominantly by 
the diffusion and therefore coincident with that of CR nuclei 
(see Figure 5 in ( [Di Bernardo et al. 2013] )). 

Our first argument aiming to constrain the value of Zh is 
the following one: for a given propagation set up, the syn¬ 
chrotron flux depends, from Eq. only on the random 
field normalization B^ani^ = 0)’ dominant respect to the reg¬ 
ular one, and on Zh, the scale-height of the diffusion region 
(X f B{r)ne{r)ds). As it is possible to appreciate in 
the Figure the fit of radio data suggests a tight relation 
-B?an(^ = 0) OC 


Secondly, we compare the observed latitude profile of the 
synchrotron emission at 408 MHz to that calculated for the 
KRA set up, setting different values for Zh- For each Zh we 
tune the value of 5ran(^ = 0) so that the average spectrum 
in these regions is reproduced (Figure]^. Low values of Zh 
are disfavoured: a analysis showed that Zh <2 kpc are 
excluded at 3cr level. 


4 Future directions and conciusions 

To fully observe and understand the GMF, there is a long way 
to go. With our analysis, we exploited the SDR as a way to 
measure the low energy LIS spectrum of CRE. For the first 
time, we have placed a constraint on the CR diffusive halo 
scale height, based on the comparison of the computed syn¬ 
chrotron emission intensity with radio observations. More¬ 
over, we stress out that - for the first time in this framework 
- our modelling of the SDR emission accounts for the pres¬ 
ence of the charge-symmetric extra-component, required 
not only to consistently model PAMELA and AMS-02 high 
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Synchrotron Spectral index Synchrotron Spectral index 



Figure 5. Spectral index map between 408 MHz and 23 GHz, with 
the assumption of a smooth 2D CR sources. Given the total inten¬ 
sity I{iy) oc (y ^ the map has been computed according the 
standard formalism: /3s = 0.248log(/ 23 // 408 ). 


energy data, but also to correctly estimate the e~ source spec¬ 
trum from CRs and radio data. In our opinion, the combina¬ 
tion of high precision CRE data, and current radio observa¬ 
tions, can be a viable method to disentangle the contribution 
of the extra-component to the total synchrotron spectrum. 

One of the greatest challenges of observing the Cosmic 
Microwave Background (CMB) in the [20 ^ 200] GHz range 
resides in the separation between the CMB and the superim¬ 
posed foreground emission: free-free, synchrotron, thermal 
dust. We pointed out that transport of charged relativistic par¬ 
ticles, and magnetic fields models should be studied simul¬ 
taneously, because both have infiuence on the synchrotron 
modelling. Synchrotron spectra may reveal signatures of spa¬ 
tially inhomogeneous particle source distributions and mag¬ 
netic fields. In our treatment here, instead of a smooth CR 
distribution invariant for rotations about the Galactic disc 
axis, we rather calculate self-consistently the synchrotron 
maps emitted by electrons whose spectral density is inho¬ 
mogeneous, due to all the relevant energy losses sustained 
while traversing regions with different distributed gas, mag¬ 
netic and radiation fields. 

Finally, the frequency range of GAMMASKY synchrotron 
simulations, from 0(10) MHz to 0(100) GHz, covers radio 
telescopes such as Planck, LOFAR, and SKA. We have been 
developed Gammas KY with the aim to support the scientific 
exploitation of the data provided by these experiments. 


References 

Di Bernardo, G., Evoli, C., Gaggero, D., Grasso, D., and Maccione, 
L.: Cosmic ray electrons, positrons and the synchrotron emis¬ 
sion of the Galaxy: consistent analysis and implications, JCAP, 
3, 036, doi: 10.1088/1475-7516/2013/03/036 2013. 

Evoli, C., Gaggero, D., Grasso, D., and Maccione, L.: Cosmic ray 
nuclei, antiprotons and gamma rays in the galaxy: a new diffusion 
model. Journal of Cosmology and Astroparticle Physics, 10,018, 
doi: 10.1088/1475-7516/2008/10/018 2008. 


Figure 6. As in Figure]^ but with the important assumption of a 
3D spiral arm structure. 


Gaggero, D., Maccione, L., Di Bernardo, G., Evoli, C., and Grasso, 
D.: Three-Dimensional Model of Cosmic-Ray Lepton Propaga¬ 
tion Reproduces Data from the Alpha Magnetic Spectrometer on 
the International Space Station, Physical Review Letters, 111, 
021102, doi:10.1103/PhysRevLett.lll.021102 2013. 

Gaggero, D., Maccione, L., Grasso, D., Di Bernardo, G., and Evoli, 
C.: PAMELA and AMS-02 e+ and e“ spectra are reproduced by 
three-dimensional cosmic-ray modeling. Physical Review D, 89, 
083007, doi:10.1103/PhysRevD.89.083007 2014. 

Ghisellini, G., Guilbert, P. W., and Svensson, R.: The synchrotron 
boiler, Astrophy. Journal Lett., 334, L5-L8, doi: 10.1086/185300 
1988. 

Gorski, K. M., Hivon, E., Banday, A. J., Wandelt, B. D., Hansen, 
F. K., Reinecke, M., and Bartelmann, M.: HEALPix: A Frame¬ 
work for High-Resolution Discretization and Fast Analysis of 
Data Distributed on the Sphere, Astrophy. Journal, 622,759-771, 
doi: 10.1086/427976 2005. 

Han, J. L., Ferriere, K., and Manchester, R. N.: The Spatial Energy 
Spectrum of Magnetic Fields in Our Galaxy, Astrophy. Journal, 
610, 820-826, doi: 10.1086/421760 2004. 

Haslam, C. G. T., Salter, C. J., Stoffel, H., and Wilson, W. E.: A 408 
MHz all-sky continuum survey. II - The atlas of contour maps, 
Astron. & Astrophys., 47, 1, 1982. 

Jansson, R. and Farrar, G. R.: A New Model of the Galactic 
Magnetic Field, Astrophy. Journal, 757, 14, doi: 10.1088/0004- 
637X/757/1/14 2012. 

Longair, M. S.: High Energy Astrophysics, 2011. 

Orlando, E. and Strong, A.: Galactic synchrotron emission with 
cosmic ray propagation models, MNRAS, 436, 2127-2142, 
doi: 10.1093/mnras/sttl718 2013. 

Pshirkov, M. S., Tinyakov, P. G., Kronberg, P. R, and Newton- 
McGee, K. J.: Deriving the Global Structure of the Galac¬ 
tic Magnetic Field from Faraday Rotation Measures of Extra- 
galactic Sources, Astrophy. Journal, 738,192, doi: 10.1088/0004- 
637X/738/2/192. 2011. 

Rybicki, G. B. and Lightman, A. R: Radiative Rrocesses in Astro¬ 
physics, 1986. 

Strong, A. W., Orlando, E., and Jaffe, T. R.: The interstellar cosmic- 
ray electron spectrum from synchrotron radiation and direct mea¬ 
surements, Astron. & Astrophys., 534, A54, doi:10.1051/0004- 
6361/201116828 2011. 






